Structural relaxation effects on interface and transport properties of Fe/MgO(001) 

tunnel junctions 
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The interface structure of Fe/MgO(I00) magnetic tunnel junctions predicted by density func- 
tional theory (DFT) depends significantly on the choice of exchange and correlation functional. 
Bader analysis reveals that structures obtained by relaxing the cell with the local spin-density ap- 
proximation (LSDA) display a different charge transfer than those relaxed with the generalized 
gradient approximation (GGA). As a consequence, the electronic transport is found to be extremely 
sensitive to the interface structure. In particular, the conductance for the LSDA-relaxed geometry 
is about one order of magnitude smaller than that of the GGA-relaxed one. The high sensitivity 
of the electronic current to the details of the interface might explain the discrepancy between the 
experimental and calculated values of magnetoresistance. 
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Tunnel magneto-resistance (TMR) is the change of the 
electric resistance of a magnetic trilayer made of two fer- 
romagnetic metallic electrodes separated by an insulat- 
ing spacer, when the mutual orientation of the electrodes' 
magnetizations changes from parallel alignment (PA) to 
antiparallel alignment (AA). The figure of merit for the 
effect is the TMR ratio, defined as 
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where Gpa (Gaa) is the conductance for the PA (AA) 
configuration. Among the possible material combina- 
tions Fe/MgO(100) magnetic tunnel junctions (MTJs) 
have attracted much attention, since extremely large 
TMR ratios were theoretically predicted on the basis of 
symmetry-driven spin filtering from the tunnel barrier—. 
These have been now demonstrated experimentally^^ 
and the symmetry filtering effect is now widely accepted. 
However, the measured TMR ratios are usually signifi- 
cantly lower than those predicted by ab initio methods. 

The origin of such a discrepancy is currently a matter 
of debate. On the theory side, interface resonance states 
(IRSs) located around the Fermi energy (Ef) are impor- 
tant for the zero-bias transport^ and are usually difficult 
to describe accurately. Moreover the MgO band-gap is 
significantly underestimated by LSDA and GGA, so that 
the barrier height might not be accurately predicted. On 
the experimental side, the quality of Fe/MgO MTJs de- 
pends on the preparation methods. The 4% lattice mis- 
match between Fe and MgO produces dislocations and 
significant relaxation at the interface. In addition, in 
typical growth conditions the Fe layers close to the inter- 
face might be partially oxidized^ 

Several different geometries for the Fe/MgO inter- 
face have been proposed theoretically, based on relax- 
ing [Fc]„/[MgO] m supercclls, consisting of n Fe and m 



MgO monolayers (MLs) subject to various constraints. 
In Ref. [U LSDA relaxation for a [Fe] 5 /[MgO] 5 superlat- 
tice, in which the in-plane lattice constant is held at the 
LSDA value for bulk Fe, yields a Fe-0 distance of 2.17 
A. In contrast an early LSDA calculation for one Fe ML 
on MgO substrate predicted an Fe-0 distance of 2.3 Ai£ 
Likewise a GGA study for one MgO ML on a Fe slab re- 
turns a Fe-0 distance of 2.21 A, and a separation between 
the first and the second Fe ML 6% smaller than that of 
bulk FefS How relevant calculations for single MLs on sur- 
faces are for Fe/MgO MTJ is however not clear and there 
is still disagreement between experiments^ and theory^. 
These controversies over the correct interface structure 
bring the question of how the different interface geome- 
tries affect the electronic and transport properties. Our 
Letter aims at answering this question. 

We perform structural relaxations for a [Fe]io/[MgO]e 
supcrlattice with both the LSDA and GGAi^ functionals 
by means of the PAW method^ as implemented in the 
VASP codei^. Charge transfer and magnetic moments 
are calculated for the differently relaxed structures with 
Bader analysis 1 ^ using the LSDA densities. A 12 x 12 
in-plane k-point grid and an energy cutoff of 400 eV 
are employed to converge the charge density. As bench- 
mark we find for bulk Fe a lattice constant of 2.75 A and 
a bulk modulus of 2.68 Mbar, both in good agreement 
with other all-electron LSDA calculations. 14 Then, the 
transport properties are computed with the SMEAGOn 1 ^ 
code, which combines the non-equilibrium Green's func- 
tions (NEGF) formalism with DFT.— The GGA PBEii 
functional and a 400 x 400 k-point grid are used for cal- 
culating the transmission coefficient. 

The layer spacings at the Fe/MgO interfaces obtained 
after structural relaxation are summarized in Tab. I. 
These have been obtained by keeping fixed the in-plane 
lattice constant, set equal to the experimental value for 
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TABLE I: (Color online) The interface structure of a [Fe]io/[MgO]e superlattice obtained by LSDA and GGA relaxation. The 
structural parameters for the unrelaxed structure are taken from Ref. [l[. The middle two Fe layers are frozen at a separation 
of 1.43 A and the experimental lattice constant for bulk Fe is used for the in-plane lattice constant. d n indicates the separation 
(in A) between layers from Fe/MgO the surface (n = 0), with the index n > (n < 0) labels MgO (Fe) layers. 



Fe, and the atomic position of the two most internal Fe 
MLs. In the table we report the obtained layer spacings 
along the stack direction d n , where do represents the Fe-0 
distance at Fe/MgO interface and n > (n < 0) refers to 
separation between MgO (Fe) MLs. In the table we also 
report the unrelaxed structural parameters taken from 
reference^. 

LSDA relaxation gives a Fe-0 distance of 2.0 A. This 
is significantly different from 2.16 A of the unrelaxed 
structure. The LSDA relaxation also predicts that the 
spacing, d—i, between the first two Fe layers next to the 
Fe/MgO interface is drastically reduced to 19% of that of 
bulk Fe. The changes in layer spacing of the MgO barrier 
on the other hand are relatively small. Interestingly, the 
GGA results are significantly different from the LSDA 
ones. In particular, GGA returns a Fe-0 distance of 2.22 
A and the reduction of d—± is less pronounced than that 
for the LSDA. Our GGA-relaxed structure is in good 
agreement with other GGA studies&i& and rather close 
to the reference unrelaxed interface structure. 
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FIG. 1: (Color online) Charge transfer, AQ, for a 
[Fe] io /[MgO] 6 superlattice as calculated from Bader analy- 
sis. Net charge transfer is obtained by subtracting the atomic 
charges of bulk Fe and bulk MgO from the Bader atomic 
charges in the Fe/MgO superlattice. The vertical dashed line 
denotes the Fe/MgO interface. 

Fig. Q] demonstrates that the charge transfer between 
Fe and MgO is very sensitive to the interface structure. 
The Bader analysis for the LSDA-rclaxed geometry pre- 
dicts that the first two Fe MLs at the Fe/MgO interface 



acquire a similar amount of electron charge, AQ, while 
for the GGA-relaxed and unrelaxed structures it indi- 
cates that the interface Fe ML loses electrons and the 
second Fe ML gains a significant amount of charge. In all 
cases MgO loses electrons to Fc. Similar charge transfer 
from MgO to neighboring metallic layers was predicted 
before for Rh/MgO interfaces^ For the LSDA-relaxed 
structure we predict a net electron transfer from MgO to 
Fc of up to about 0.06 electrons per atom. This is larger 
than the one for the other two geometries. The Bader 
charges for O atoms at the interface and in the middle 
of the MgO barrier are respectively 7.64 and 7.71. This 
indicates that charge is transferred to Fe mainly from O. 
Finally, since the GGA-relaxed structure is close to the 
unrelaxed one, the charge transfer is also similar. 

The local magnetic moments of the Fe atoms [Fig.[2Ja)] 
are calculated by integrating the spin densities in the do- 
mains determined by charge densities resulting from the 
Bader analysis. These are similar for the GGA-relaxed 
and unrelaxed structures, both presenting a significant 
enhancement of the Fe magnetic moment at the MgO 
interface. For the LSDA-relaxed structure however such 
interfacial magnetic moment is dramatically suppressed. 
This is only 1.10 /is, to be compared with 2.65 \ib and 
2.68 [Ib respectively for the GGA-relaxed and unrelaxed 
geometries. Since AQ is much smaller than the change 
in the magnetic moment, the dramatic decrease of inter- 
face magnetic moment in the LSDA calculations is caused 
by an internal electron redistribution between the major- 
ity and minority spin sub-bands. This is demonstrated 
by the spin-decomposed electron occupation [Fig. |2^b)]. 
Since in the LSDA relaxation <i_i is significantly shorter 
than that for GGA relaxation, the suppression of the 
magnetic moment is expected, based on the fact that 
many magnetic materials undergo magnetic collapse un- 
der pressured 

The effects of the interface structure on the electronic 
transport are studied with SMEAGOLi^ 5 . At zero bias the 
conductance, G, is simply G = e /hT(Ep), where e is 
the electron charge, h the Planck's constant and T(E-p) 
the transmission coefficient calculated at E-p. T(Ep) is 
presented in Fig. [3] The transmission for majority spins 
in the PA depends weakly on the energy but it is re- 
duced by about one order of magnitude for the LSDA- 
relaxed structure, when compared to both the GGA and 
the unrelaxed ones. In contrast the minority spin chan- 
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FIG. 2: (Color online) The magnetic moment (a) and spin- 
decomposed electron occupation (b) for a [Fe]io/[MgO]e su- 
perlattice. These are calculated using Bader analysis for the 
valence electrons. The vertical dashed line in the top panel 
denotes the Fe/MgO interface. 
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FIG. 3: (Color on line) Transmission T(E) for PA majority 
spins (a), PA minority spins (b), and for the A A (c) of a 
Fe/MgO/Fe MTJ. The Fermi energy is set to zero. 



nel for PA is dominated by a peak at around Ep. This 
originates from IRSs at the Fe/MgO interface.— For the 
LSD A- relaxed structure these peaks are shifted 0.2 cV to 
lower energies with respect to the GGA-relaxed case. A 
similar shift is found for both spins in the AA configu- 
ration. This shift leads to a large change in the low bias 
conductance. As is shown in Tab. [Til the total conduc- 
tance of the LSDA-relaxed structure is about 30 times 
smaller than that of GGA-relaxed structure for the PA. 
Although the magnetic moment and charge transfer are 
similar for the GGA-relaxed structure and the unrelaxcd 



structure, the transmission coefficients for the two struc- 
tures are significantly different. The GGA-relaxed struc- 
ture shows less pronounced and wider peaks in the trans- 
mission coefficient than those of the unrelaxed structure. 
This results in the conductance of GGA-relaxed structure 
being 5 times smaller than that of the unrelaxed one. 
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FIG. 4: (Color on line) The projected density of states of the 
interface Fe 3d orbitals in a [Fe]io/[MgO]e junction. Posi- 
tive and negative projected DOS represent the majority and 
minority components, respectively. 

Since the peak in transmission at around Ep originates 
from the IRSs at the Fe/MgO interface, the shift in en- 
ergy of the transmission coefficient for the LSDA-relaxed 
structure relative to that of the GGA-relaxed one is the 
result of a high sensitivity of IRSs to the interface geom- 
etry. In Fig. [4] the density of states (DOS), projected on 
the interface Fe layer, is shown for both spins. For the 
LSDA-relaxed structure this is shifted by about 0.2 eV 
with respect to that of the GGA-relaxed structure, which 
causes the corresponding energy shift of T(E). Since 
the unrelaxed and GGA-relaxed structures are close to 
each other, the general features of the projected DOS 
for the two are also similar. In order to get a better in- 
sight into the role of relaxation, in Fig. [5] we show the k- 
resolved projected DOS at E-p for the interfacial Fe layer. 
The results for the unrelaxed and GGA-relaxed struc- 
tures are rather similar, but they differ remarkably from 
that for the LSDA-relaxed case. In general the LSDA- 
relaxed geometry shows a substantial reduction of the 
spin-polarization, demonstrated by the fact that the k- 
resolved DOS is similar for the majority and minority 
spins. This is expected to produce a reduction in spin- 
filtering and as a consequence a reduction in TMR. It is 
also important to note that the main differences between 
LSDA-relaxed and GGA-relaxed (and unrelaxed) geome- 
tries are more evident around Brillouin zone center, i.e. 
for electrons with a large linear momentum perpendicular 
to the MgO barrier and therefore with a larger transmis- 
sion probability. 

The changes in conductance lead to very different val- 
ues for the TMR. Our calculated values for the TMR 
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FIG. 5: (Color on line) The k-resolved projected density 
of states for the interface Fe ML in a Fe/MgO/Fe junction. 
Unrelaxed, LSDA-relaxed and GGA-relaxed structures are 
shown in the top, middle and bottom panels respectively. The 
left and right panels are for majority and minority spins, re- 
spectively. 



TABLE II: The TMR and zero bias conductance for PA and 
AA a Fe/MgO/Fe MTJ. GGA PBE is used for all transport 
calculations. The results are shown for the unrelaxed, the 
LSDA-relaxed and for the GGA-relaxed structures. The ex- 
perimental TMR is for a Fe/MgO/Fe MTJ with a 2.3 nm 
thick MgO barrier and a temperature of 20 K.— 
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are shown in Tab. [TTJ The LSDA-relaxed structure has 
a TMR of about 10 3 , which is much smaller than the 
TMRs of the other two structures, and in much closer 
agreement to experiments^ The dramatic reduction of 
TMR for the LDA-relaxed structure originates from two 
features: the shift to lower energies of the IRSs trans- 
mission peak, and the large reduction of the majority 
spin transmission for PA. These can be both associated 
with a loss of spin-polarization of the first Fc layer at the 
Fc/MgO interface in the LSDA-relaxed case. 



In conclusions, motivated by the significant difference 
between theory and experiments on the reported TMR 
values in Fe/MgO/Fe(100) MTJs, we studied the inter- 
face structure and its effects on the electronic and trans- 
port properties by means of DFT and the NEGF for- 
malism. In general LSDA- and GGA-relaxed interfacial 
geometries are rather different, yielding different charge 
transfer and interfacial magnetism. In particular the lo- 
cal magnetic moment of the interfacial Fe ML is severely 
suppressed for the LSDA-relaxed structure and enhanced 
for GGA-relaxed ones. These differences are reflected in 
the transport properties. In particular the differences at 
the interface determine the energy position of resonances 
in the transmission for the minority spins. These reso- 
nances are caused by IRSs and largely determine the zero- 
bias conductance. As a consequence the TMR is rather 
sensitive to the interfacial structure. These features can 
partially explain the disagreement between theory and 
experiments and are expected to apply to many other 
systems relevant for spintronics 
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